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ABSTRACT 



Monterey Bay, on the west coast of the United States, is unique in 
that it is a large, symmetric, semi-elliptical bay divided north and 
south by the deep Monterey Canyon. The effect of the canyon on seiching 
within the bay and on long wave oscillations within the bay was studied 
by analyzing synchronized wave records at each end of the bay. Power 
spectra and cross spectra calculated for five periods selected from six 
months continuous data indicate the Monterey Canyon has a profound effect 
on the bay’s oscillating characteristics. The canyon appears to act as 
an impedence barrier dividing the bay into two independent oscillating 
basins each having recurring long-period waves which persist during 
significant long-wave activity. 
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I. INTRODUCTION 



Monterey Bay (Figure I) lies on the west coast of the United States, 
located about sixty nautical miles south of San Francisco. The Bay is 
a large semi-elliptical bay which has some topographical features which 
are unique. A deep trough of the North Pacific basin, bounded on the 
north by the Mendocino Seascarp and on the south by the Murray Sea- 
scarp, approaches closer to the coast at Monterey Bay than at any other 
point along the North American coastline. As a result, long waves propa- 
gated across the Pacific Ocean in an easterly direction may be somewhat 
contained. The continental shelf within the bay extends out to approxi- 
mately the 600 foot depth contour and is cut by the deep Monterey Canyon 
which has a volume of over 50 cubic miles. Monterey Canyon will be seen 
to exert significant influences on the oscillating characteristics of 
the Bay. 

There have been several studies concerned with the effects of long- 
period waves in Monterey Bay (a long-period wave for this study is 
defined as a wave having a period in excess of 1 minute), but these 
studies have not been able to define the causes of the long wave 
activity and have left many areas of inquiry. Hudson [1947] proposed 
a surge model of Monterey Harbor for the U, S. Army Corps of Engineers. 

He used six months data (October 1946 - April 1947) from three elec- 
tronically synchronized sea-surface recorders. These results are 
tabulated in Table I. The presence of the long-period waves was 
observed. However, the cause of the surge mechanism was unknown. 

Wilson [1965] conducted a long -wave study of Monterey Bay for 
the U. S. Army Corps of Engineers to determine the feasibility of 
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construction of an engineering model of the surge phenomenon that 
occurred at various times within Monterey Bay. Statistical data were 
collected from sensors located at Monterey Harbor, and Santa Cruz 
(Figure 1). Three long-period wave recorders were also installed 
within Monterey Harbor and were in continuous operation from October 
1963 to April 1964. Two were arranged so that tides and sea-swell 
were filtered out, and a third recorded swell approaching the harbor. 
Wilson concluded that the surge within Monterey Harbor was not a result 
of incoming swell as no correlation was found to exist between long- 
period waves within Monterey Harbor and the incoming sea-swell. The 
Santa Cruz sensor was designed so that it functioned as a long-wave 
recorder. The recorder operated continuously from October 1963 to 
February 1964. Wilson’s results for Monterey Harbor are summarized 
in Tables II and III. 

Wilson performed a residuation analyses of different records from 
sensors located around the bay in order to determine a local evaluation 
of the oscillations within the bay. Residuation analyses is accomplished 
by successively subtracting "apparent" periods from the wave record 
until a smooth trace remains. The subjectivity of which this method is 
accomplished can lead to ficticious results. These results are listed 
in Table IV. Table V is a synopsis of spectrum analyses for three days 
record of the Monterey sensors. Wilson concluded, based on Tables IV 
and V along with calculated modes of oscillation, Figures 24-29, that 
Monterey Canyon functions as an impedence between north and south 
portions of the bay, and therefore, free oscillations are to a large 
extent uncoupled. Wilson further states that the effect of any sharp 
discontinuity in submarine topography, equivalent to the edge of a 
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continental shelf, is to serve as a nodal position for any shelf 
oscillation to which the shelf is susceptible. 

Robinson [1969] analyzed long-wave activity in Monterey Bay on 23 
January and 20 April 1969. Tide recorders were located at Monterey 
Wharf #2 and Santa Cruz Municipal Wharf, however, the records were not 
synchronized. Thus, cross spectra and phase information was not 
obtainable. Robinson concluded that Monterey Canyon affected the 
seiching motions within the bay, however, it did not appear to divide 
the bay into two independent oscillating basins, as concluded by Wilson 
[1965]. Similar periods and amplitudes were found at Monterey and 
Santa Cruz concurrently, but correlation between the two locations 
could not be computed due to non-synchronization of records . 

In the present study simultaneous tidal records from tide gages 
located at Monterey Wharf #2 and Santa Cruz Municipal Wharf were 
analyzed during periods of significant long-wave activity in Monterey 
Bay. Continuous analog records were maintained from 3 November 1969 
to 30 April 1970 and five days (6 November, 3 December, 6 December, 

14 January, 20 January) were selected for analysis and comparison. 

The objectives were to perform an overall long-wave study of 
Monterey Bay, and to examine various modes of free oscillation of the 
bay to determine if in fact seiching does exist in the bay. This was 
accomplished through examination of the individual power spectra at 
Monterey and Santa Cruz, computation of the cross spectra and deter- 
mination of phase differences of long-period waves present. 
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TABLE I 



Long-Period Waves in Monterey Harbor (after Hudson, 1949) 



PERIOD 


AVERAGE HEIGHT 


PERCENTAGE OF 


(min) 


(ft) 


TIME PRESENT 


1-2 


0.4 


20 


2-4 


0.5 


30 


4-15 


not given 


15 



TABLE II 

Marine Advisor's Data for Monterey Harbor (after Wilson, 1965) 



SENSOR 


PERIOD 


AVERAGE HEIGHT 


PERCENT OF 






(min) 


(ft) 


TIME PRESENT 




1 


1.7-14 


0.1-2. 5 


0-50 




2 


1-14 


0. 1-3.0 


0-55 






* 


TABLE III 




Marine Advisor 


f s Data 


for Santa Cruz Harbor 


(after Wilson, 1965) 


PERIOD 




AVERAGE HEIGHT 


PERCENT OF 


(min) 




(ft) 


TIME PRESENT 


1-2.3 


n 




0-18 


2.3-4 
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| 

^ 0. 2-2.0 


20-70 


10-14 
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0-25 
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Periods of Oscillation from Resi dilation Analysis (From Wilson, 1965) 
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Periods of Oscillation from Spectral Analysis (From Wilson, 1965) 



II. INSTRUMENTATION 



I 

A. kONTEREY 

The data for Monterey was obtained using a standard Coast and 
Geodetic Survey automatic tide gage. (Manual of Tide Observations, 
1965), located on Monterey Municipal Wharf #2, Monterey Harbor. The 
tide gage is maintained daily by NPS personnel. 

The Monterey tide gage senses changes in water level by means of 
a float/pulley arrangement, which is entirely mechanical. The recording 
drum is advanced by a clock mechanism. The drum speed is designed for 
1 in/hr, however, the hourly feed was measured to be 1.03 in/hr. The 
instrument is time checked daily with time marks being accurate within 
+ 5 seconds. The marigram is recorded in rectilinear coordinates on 
plain white paper. A stilling well, which is a 12 inch diameter steel 
pipe with a 1 inch orfice in the bottom, serves as low pass filter for 
the majority of wind waves (periods of one minute and below). 

B. SANTA CRUZ 

The Santa Cruz data were obtained with a Bristol Model 28 gas- 
purging pressure (bubbler), portable tide gage, (Figure 2), located 
on Santa Cruz Municipal Wharf. This instrument senses changes in 
water level by means of a nitrogen-filled tube which is connected to 
a bellows system (Manual of Tide Observation, 1965). 

A bubbler orfice chamber was connected to the end of the sensing 
tube to reduce sensitivity to short period wave action. There is also 
a bellows inlet needle valve which could be throttled to filter add- 
itional wave action from the record. 
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1. NITROGEN BOTTLE 

2. PRESSURE REDUCING VALVE 

3. PRESSURE DIFFERENTIAL REGULATOR 

4. ADJUSTABLE NEEDLE VALVE 

5. TRANSPARENT BUBBLER CHAMBER 

6. PULSATION DAMPER 

7. STRIP CHART RECORDER WITH TRANSDUCER 

8. TELEMETERING PRESSURE TRANSMITTER 

9. TUBING 

10. BUBBLER ORIFICE CHAMBER 

11. STILLING WELL 



FIGURE 2 THE BUBBLER TIDE GAGE. 





The pressure reduction valve is regulated to allow a constant 
pressure greater than that due to the maximum head anticipated over 
the orf ice . The pressure differential regulator together with the 
shunt line provide for a constant pressure difference of about 3 
pounds per square inch across the flow regulator. Thus, the rate of 
flow can be adjusted to a desired constant value and is not dependent 
on the tide stage. Flow is through a silicone oil-filled bowl so that 
it can be monitored visually. A rate of 30-60 bubbles per minute was 
used for this study. 

Advantages of this scheme are that the sensitive elements of the 
gage need not be designed to withstand the underwater environment; and 
fluctuations in barometric pressure are not directly reflected on the 
record . 

A significant portion of the system capacity is at the submerged 
orf ice , thus providing some high frequency filtering and protecting 
the small diameter supply line and orf ice from marine fouling. 

A Bristol Company clock recorder provides an analog record. The 
record is on a 6 inch (15.2 cm) wide strip chart and is recorded in 
curvilinear coordinates. Paper advance is adjustable and 6 in/hr was 
utilized. The clock drive has an eight-day spring, and the record is 
sufficient for about seven days. Manufacturer's claims are that 
hysteresis and/or non-linearity are limited to Vjo of the full-scale 
value. For sea water use, the instrument is calibrated to correspond 
to a specific gravity of 1.025. 

In addition to the above high-frequency filtering mechanisms, a 
stilling well was designed and installed to attain the desired 
filtering of high frequency and waves [Robinson, 1966]. The well was 
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constructed of a 20 foot section of polyvinyl chloride (pvc) 6 inch 
inside diameter pipe. It was capped on 16, 1/4 inch inside diameter 
holes drilled in the side. Copper sleeves were inserted to eliminate 
fouling. The 16 holes provided the capability of increasing the orfice 
from a 1/4 inch to a 1 inch diameter opening. 

The response characteristics of the well were determined theoretic- 
ally for two orfice sizes and three different wave frequencies using 
the equation for the rate of water rise in a well [Robinson, 1969] . 

dh^/dt = 0.6 a/A - h^) 

where , 

a = orfice area 
A = well area 

0.6 = empirical orfice flow coefficient 

g = acceleration due to gravity 

hQ = water height outside the well, i.e., the forcing function 

h^ = water height insdie the well. 

The forcing function, h^, was chosen to be a simple sine function 

of unit amplitude, and of frequency equal to the wave frequency of 

interest. The initial conditions were h. = 0 at t = 0. The results 

i 

are summarized in Table VI. 

The response characteristics for the Monterey stilling well were 
not calculated since the orfice area to well ratio is larger than the 
Santa Cruz well, providing response characteristics better than those 
of the Santa Cruz. A 1/4 inch orfice was used in the Santa Cruz well. 
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TABLE VI 

Computed Response Characteristics of Santa Cruz Stilling Well 

(from Robinson, 1965) 



PERIOD 


ORFICE DIAMETER 
(in.) 


PHASE LAG 
(deg) 


AMPLITUDE 

REDUCTION 

(%) 


20 sec 


0.25 


180 


95 


20 sec 


1.00 


72 


45 


60 sec 


0.25 


75 


91 


60 sec 


* 

1.00 


30 


5 


25 min 


0.25 


6 


1 


25 min 


1.00 


0 


0 
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III. ANALYSIS PROCEDURE 



A . MONTEREY 

The Monterey tidal records were digitized at 100 points per inch 
giving a sampling interval of 34.95 seconds. The discrete Fast Fourier 



analyzed, 2048 data points were used which gave a record length of 
71545 seconds, or approximately 20 hours. In order to remove the tidal 
influence and reduce "leakage", the records were high pass filtered by 
fitting a least-squares polynomial curve to the raw data and the curve 
was then subtracted from each ordinate point (Appendix C) . An example 
of the fit to sample data points is shown in Figure 3. It was found 
that the 6th-degree polynomial provided the best fit for this particular 
tidal cycle. 

Once the tidal effects were removed, the data is subjected to the 
Fast Fourier Transform to obtain the energy density and phase estimates 
as a function of frequency. The IBM/360 standard subroutines HARM and 
RHARM were used to calculate the Fourier coefficients, from which the 
one -dimensional spectra are derived. This method and procedure is 
explained in Appendix B. 

Subroutine RHARM gives the raw Fourier’s coefficients 



Transform utilized requires 2x2 m data points, and for the records 




which are then combined as 




A 



iB 



2 



2 



j = 1,2, ...N 



to give the one-sided spectral estimate 



23 




seconds ) 




Figure 3 

Least Squares Curve for Santa Cruz and Monterey, 6 Dec 69 
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where 



v (V-/)© 



T = record length (seconds) 

Ej = spectral estimate (m^-sec) 
The associated frequencies are 




The raw spectral estimates were then smoothed over five frequencies, 
giving each new estimate 10-degrees of freedom and increased confidence 
in the spectral peaks. Smoothing over five spectral estimates removes 
spurious peaks, but does not smear the actual peaks. The 90% con- 
fidence limits for 10 degrees of freedom are interpreted as meaning 
that one can state with 90$ accuracy that the actual spectral estimate 
is greater than . 49E^ but less than 1.60Ej, as the power spectrum is 
distributed according to a CHI-SQUARE distribution. 

B . SANTA CRUZ 

The analysis procedure of the Santa Cruz data was more complex than 
the Monterey data due to two considerations : 

(1) The Santa Cruz tidal trace is recorded on curvilinear 
coordinates, but is accepted as rectilinear coordinates 
by the CAIMA 480 digitizer. 

(2) The sampling interval of Santa Cruz is 6 seconds which 
must be interpolated and matched point for point in time 
with the Monterey data, in order to compute the cross 
spectra . 
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The unadjusted amplitudes and time increments were obtained from use 
of the CALMA 480 digitizer, and computer program CONVERT. The analog 
recoird was recorded on Bristol Company chart 6112, shown in Figure 4, 
which has curvilinear coordinates (t,R0). The digitizer then recorded 
incremental pairs, in the order X and Y (of the cartesian coordinates 
of the tidal curves), every .01 inch of stylus travel. It was there- 
fore necessary to form the cartesian coordinates by summation within 
the computer. These cartesian coordinates were then converted to the 
approximately correct T, R0 coordinates by employing geometric 
relationships . 

The computation of the correct value of t^ (time axis) is critical 
to the conversion. On the time scale, 0.1 inch of record was equivalent 
to a sample interval of 6-seconds. The summation of X, X^ , representing 
the horizontal travel was converted to the t, R0 coordinates by 




l 



l 



l 



where , 



t. = time in curvilinear coordinates 

l 



R = radius of curvature for arcs of constant t. 



i 



Y. = the ordinate on the cartesian scale 



l 



X^ = time scale in equally spaced increments of t on the 



cartesian scale. 



These geometric relations are shown in Figure 5 where 




The adjusted heights are then computed as: 



Y. = R0 



l 



where Y^ is now the adjusted sea-surface elevation value. Since the 
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Bristol Company Chart 6112 
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Geometric Relations for Conversion from Curvilinear Coordinates (R,9) 

to Rectilinear Coordinates (Y,t) 



space increment is .01 inch, the adjusted time, t. is 



t . 
1 



= 600 X. 

i 



The adjusted sea surface elevation (with reference to an arbitrary 
level) and adjusted times, were not in general equally spaced. There- 
fore, an interpolation scheme was used to obtain data points at a 
sampling interval of 34.95 seconds in order to match the Monterey data, 
point for point (Figure 6) . 

A linear interpolation was performed between time increments greater 
than 34.95 seconds. The slope of the line is 



the desired interval of 34.95 seconds is exceeded, then interpolates 
between the value greater than 34.95 seconds and the prior value which 
has an increment less than 34.95 seconds. Santa Cruz data were then 
aligned with the Monterey data point for point. 

The first 2048 points of the Santa Cruz data were handled in the 
same manner as the Monterey data. The tidal influences were removed 
by the least-squares scheme, and the adjusted data were analyzed in 
the same manner as the Monterey data. 




The horizontal spatial scale for the interpolated values is, 



X. = (i-1) A 

l 



where A = 34.95 sec. 



The resultant interpolated amplitude is, 




Effectively, the computer checks each successive time value until 
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Figure 6 

Linear Interpolation 



C. CROSS SPECTRA 

Given two random, stationary time-series, f^(t) and f ^ (t) the 
cross-spectra is defined as: 

CD 

$ 12 (o) = hi { ^12 ( T ) e " 1GTdT (!) 

-CD 

where , 

T 25 time lag 
O = angular frequency 

cp^2 ~ the cross-covariance between f^(t) and f^Ct) gi ven by 

T/2 

cp i2 < T > = lim f ^(t) f 2 (t+T ) dt (2) 

too j 

-T/2 
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where 



T = time lag 
T = .record length 
t = time 

If the time lag t is replaced by -t and substitution of the dummy 
variable 6 = t - t we have, 

T/2 

cp 12 (“T> = lim f f x (6 + t) f 2 (6)d6 
-T/2 

then, the cross -covariance is an even function such that 

Cpi2 <-t) = Cp 2 i (t) (3) 

The cross-spectrum (1) can be written (using Euler's formula) 
e " laT = COS oT-i SIN ot 



as 



*12 



a 2n\ f CP 12 ( T ) C0S OTdT-i J CP 12 ^ T ) SIN OTdTj 



since not > g en ^ral, an even function, therefore, the cross 
spectrum is a complex function of a. 

The cross spectrum combines both the cosine and sine transforms of 
cp^CT). From this, if the cross-spectrum can be resolved into odd and 
even complnents, it can then be written as, 



-h\ J" cp even (t) COS crrdT-i j cp odd (t) SIN OTdT> 



(4) 



since the cosine (sine) transform of an odd (even) function is zero. 
For convenience we define two functions A(t) and B(t) where 
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A 12 (T) = 



B 1 2 (T) = 



CP 12 (T) - Cp 21 (T) 



cp 12 (t) + cp^CO 



( 5 ) 



If T is replaced by -T in (5), and using (3) gives = ^^(t) 

and w * iere ^ an ° dd function and B even. Since 

cpf 2 (t) is resolved as 



cp 12 (t ) = a 12 (t) + B 1 2 ^ T ^ 

we see A and B are the odd and even parts respectively of cp^* 

Using this equation, (4) can be written as 

Too 00 ^ 

$ 12 (o) )f b 12 (t) cos aTdT_i / A 1 2 (T) sin 0TdT r * (6) 

t -00 -OO J 

The co-spectrum is defined as the real part of (6) 



c 12 (o) f B 1 2 (T) sin aTdT 



while the negative of this imaginary part defines the quadrature- 
spectrum 



Q 12 (o) = ■ Itt / a 12 (t) SIN aTdT * 



The cross spectrum is then written in complex form 
$ 12 (o) = c 12 (o) - iQ 12 (a) 



(7) 



If a time series is correlated with itself the resulting co- 
variance function is even. The cosine transform of the auto-covariance 
function is the power spectrum of the record. The co-spectrum is the 
cosine transform of the even part of the cross-covariance function, 
and gives the amount of power due to in-phase (0°) or out of phase + 
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180° components between the two records. The quadrature spectra is 
similar, however, the components are out of phase by a constant amount. 
The spectral phase difference defined by 



e 12 (o) = TAM 



Ql2 (0> 

C 12 (a) 



and is the argument of the complex conjugate of the spectrum. 

In this study the Fast Fourier Transform method was used to determine 
the power spectra vice the auto-correlation method. This resulted in a 
slightly different scheme of calculating the cross-spectra, as described 
below. 

From subroutine RHARM the Fourier coefficients 



^ ti 

2 ’ B 0 °’ 2 ’ 2 ’ ' " 2 ’ B N 0 






were determined for both the Monterey and Santa Cruz data. 

From the Fourier coefficients we form the modules for each station, 

F M ( j> = 2^ + ~T* I F mI G 1 (Monterey) 

and 



A iB 

F s (j) = ^ + ~~ 2 ~ ~ |Fg|e le 2 (Santa Cruz), 



or 



Fg ( j ) = F s (j)e' le 2 




(complex conjugate) 



where j = 0,1,2, ...N. 

Then the cross-correlation function, 

2N-1 . . _ 

cp^^) - S F g (j) F M (j)e 1J0 l 

j=l 

where = 2 tt/N dt 



t 
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and the cross spectra 



T/2 

*12 0)- -i / ^(Oe-^VdT 

-T/2 

form a Fourier integral pair. 

The cross-spectra $^0) can be expanded as 

* 12 0) - F M (j) Mi) - |F M (j)|e ie l|F Q (j)|e‘ iG 2 



'M' J ' S 



M' 



= |F M (j)||F Q (j)|e i(e r G 2 ) 



where 



; i - tan ' 1 0 



(Monterey) 



- tan ' 1 ( 49 - 



(Santa Cruz) 



By use of Euler f s formula, 

|A|e^ = | A | (cos 0 + i sin 0) 
we have, expanding 

* 12 (j) = |F M (j)l|F s (j)| COS ( £l - e 2 ) 

+ i|F M (j)||F s (j)| SIN ( £l - e 2 ), 



or, $^2 = [co-spectra (in phase)] + i [quad-spectra (out of phase)]. 

For the co-spectra, if COS (0) = 1, and we have maximum 

value for in phase. If = + t ^ ie < l ua ^ spectra, SIN tt/2 = 1 

and we have maximum out of phase. If we arbitrarily let = 0 then 

sin c , 

1 QUAD 

tan g . = = 

1 cos Gl CO 

and the phase difference, 

« - er> 
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From this, 
difference 



the co-spectra, quad-spectra, cross-spectra and phase 
was computed for each spectral component. 
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IV. ANALYSIS RESULTS 



I 

A. -IDENTIFICATION OF LONG WAVES PRESENT 

Spectral analysis of ten tidal records indicated that the long-wave 
activity in Monterey Bay is characterized by a few specific spectral 
peaks which dominate the spectra. Their period identification and 
possible causes are explained below. 

51.8 Minute Period 

According to Wilson (1965), the fundamental mode of longitudinal 
oscillation of the north-south extremities of Monterey Bay is 44.2 
minutes. If the bay is approximated as an enclosed basin of uniform 
depth and length, (Figure 7), the periods of oscillation are roughly 
calculated in the following manner 



Monterey 




Santa Cruz 



Figure 7 

Fundamental Longitudinal Oscillation of Monterey Bay 
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The shallow-water wave celerity is approximated by. 



c = Vgh 

where , 

C = wave celerity 
g = acceleration due to gravity 
h = average depth of bay 
The wave length of the fundamental period is 

L = 2w 

letting C = L/T, and substituting gives, 

2w 

T 1 'Vih 

or in general 

1 2w 

T n = nTgh n = 1, 2, ... 

for the harmonics of the fundamental mode of oscillation. 

Substituting , 

w = 20nm 
h = 240 feet 
g = 32 .2 ft/sec^ 

gives the fundamental period of oscillation 

T ^ = 46,5 minutes 

The 51.8 minute long-wave appears throughout the different spectra 
and is evaluated as the fundamental period of longitudinal oscillation 
36,1 Minute Period 

The fundamental mode of transverse oscillation of the east-west 
extremities of Monterey Bay is computed by Wilson [1965] to be 32.3 
minutes. Again, as in all Wilson’s calculations, he assumed a nodal 
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line from Pt . Santa Cruz to Pt . Pinos. From Figure 8 the fundamental 
transverse oscillation is calculated. 



Assumed Node 




Moss Landing 



Figure 8 

Fundamental Transverse Oscillation of Monterey Bay 
The fundamental wave lerfgth is 

L = 4b 

Using the shallow water wave equation 

C = Vgh = •£ 

and substituting gives, 

1 4b 

T „ ' nTih « - 1, 2, ... 

For Monterey Bay, 
b = 8nm 
h = 240 ft 

the fundamental period is calculated as 
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= 36.3 minutes . 



The calculated fundamental transverse mode, which is in agreement 
with’ Wilson is present in Monterey four out of the five days analyzed, 
however, it only appears in Santa Cruz on 3 December 1969. 

In addition to the above, one can also consider the possible pre- 
sence of significant edge -wave activity [Munk, Snodgrass and Carrier, 
1956]. The edge-wave is a wave form in shallow water on sloping 
shelves and beaches. It is called an edge-wave because the waves 
approach with the crest perpendicular to the beach and travel along 
the coast. The waves are sinusoidal, and it is assumed the shelf 
slope is constant. The celerity of the edge -wave as a function of 
period and shelf slope is given by, 



c - §“ Sin 3 

where , 

C = wave celerity 
T = period of edge -wave 
3 = shelf slope 

g = acceleration due to gravity. 

Substituting C = L/T we have the fundamental period and harmonics 
of the edge -wave 



D = /— -r“ (2n- 1) 
n Vg sm 3 



n = 1, 2, 



where L is the wave length equal to twice the longitudinal distance 
across the bay for the fundamental mode. 

Substituting , 

g = 32 . 2 ft/sec 
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3 = 600 ft/8nm 
~ = 20nm 

the fundamental period of oscillation and associated harmonics of edge- 
waves for Monterey Bay are computed as , 

T = 35.8, 25.3, 20.2, 17.8, 15.7, ... minutes 

These modes are also called "trapped modes" from wave-guide theory 
as the wave energy is trapped into a continental wave guide where in 
theory the energy remains trapped indefinitely under idealized conditions. 

Although the 35.8 minute long-wave is significantly present as 
indicated above, it is not evaluated as an edge-wave as phase computations 
do not support any pronounced edge wave activity. 

27.7 Minute Period 

The 27.7 minute period appears to be the most persistent long-wave 
in the bay. It appears in the spectrum having significant energy 
density on four out of five days analyzed, and is present at both Santa 
Cruz and Monterey on 6 December 1969 and 14 January 1970. The 27.7 
minute wave is also the dominant long-wave in Robinson’s [1965] analyses. 
Such waves can be interpreted as shelf-waves [Munk, 1962]. The formula 
for the case of normal incidence on a two-step topography (Figure 9) 
consisting of the continental shelf of width A and constant slope S f , 
and the continental slope of width B with slope S is 
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where , 

Tq = fundamental period of the shelf wave 
(3 = constant which is dependent upon ratios S'/S and h'/h 

where h'/h is the ratio of depth of break in continental 
shelf to depth of sea floor at the bottom of continental 
slope. 

g = acceleration due to gravity. 

For the case of Monterey Bay, choosing the 100 fathom curve as h' 
and average distance A equal to 8.0nm and letting 4g / tt = 1.54 when 
S' < < S and h' < < h, the formula for the shelf-wave fundamental 
period becomes 

x 0 (60)1.54 Jgh 7 * 

Substituting , 

A = 48000 ft (average distance to 100 fathom curve) 
g = 32 .2 ft/ sec^ 
h' = 600 ft 
we have , 

Tq = 29.6 minutes . 

Modification of A and h' in the shelf-wave formula for different 
locations around the bay changes the fundamental period only slightly. 

The remaining long -waves which are present during significant long- 
wave activity are those of periods 22.5, 18.9, 16.3, 13.6, 10.1, and 
9.3 minutes. These are considered to be harmonics of the basic funda- 
mental period 51.8, 36.1 and 27.7 minutes, and are in general weak in 
spectral density. The significant spectral peaks found on various dates 
are summarized in Table VII. 

41 
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Summary of Significant Spectral Peaks in Monterey Bay 



B. INTERPRETATION OF INDIVIDUAL SPECTRA 



Figures 10, 12, 14, 16, and 18 show the synchronized tidal traces 
of Santa Cruz and Monterey on dates indicated. The traces are plotted 
with a sampling interval of one minute which is approximately equal to 
the Nyquist period of 70 seconds, which is the lowest period which can 
be analyzed using the 34.495 sampling rate. The tidal traces are repro- 
duced here for the first 16 hours of each period analyzed. 

6 November 1969 

The individual spectra for 6 November 1969 are plotted in Figure 11. 
The record tends to be a "noisy" record with considerable long-wave 
activity. The weather conditions in Monterey Bay were calm winds and 
seas in the early morning, increasing to 15 kts of wind from the south- 
west about 2000. The surface air temperature was 55°F and skies were 
overcast for most of the day. The 42.5 minute non-recurring long 
wave appears to be significant (the term significant peak used in this 
study is a peak that has significant spectral density when compared 
with other peaks in the same record vice having a specific value) at 
both Monterey and Santa Cruz only on this date. Most noticable in the 
records are waves of periods 51.8, 42.5, 16.3 and 10.1 minutes which 
appear at both stations and have significant energy-density. 

3 December 1969 

Figure 13 shows the individual power spectra for Monterey and Santa 
Cruz from 0130-2130 3 December 1969. Examination of Figure 12 shows 
that significant wave activity in the bay commenced about 0950. 

Initially the surface winds were from the north at 2 kts, shifting to 
south at 15 kts by 1200. Skies were clear and surface air temperature 
increased from 45°F in the early morning to 58°F in the afternoon. The 
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Sea Surface Height (feet) Sea Surface Height (feet) 





Figure 10 

Sea Surface Heights, Santa Cruz/Monterey 
0730 6 November 1969 to 0330 7 November 1969 
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Frequency (millihertz) 

Santa Cruz 0730 - 0330, 6 Nov 69-7 Nov 69 




Frequency (millihertz ) 

Monterey 0730 - 0330, 6 Nov 69-7 Nov 69 



Figure 11 

Spectral Wave Analysis, 6 Nov 69 
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Figure 12 

Sea Surface Heights, Santa Cruz/Monterey 
0130 - 2130, 3 Dec 69 
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Energy Density (meters -sec) Energy Density (meters -sec) 




Frequency (millihertz) 

Santa Cruz, 0130 - 2130, 3 Dec 69 




Frequency (mi 1 liher tz ) 
Monterey, 0130 - 2130, 3 Dec 69 

Figure 13 

Spectral Wave Analysis, 3 Dec 69 
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spectra shows intense long-wave activity at Santa Cruz and moderate 
activity at Monterey. The effect of persistant southerly winds may 
account for this difference, however, no correlation is drawn. The 
36.1, 13.6 and 10.1 minute waves have pronounced energy at both 
locations, while the 51.8, 27.7 and 16.3 minute waves are significant 
at Santa Cruz alone. Note, that at Santa Cruz the 13.6 minute wave 
can very well be a harmonic of the 27.7 minute shelf wave while the 
16.3 and 10.1 minute waves could be a harmonics of the fundamental 
transverse mode of 36.2 minutes. 

6 December 1969 

The tidal traces for 6 December 1969, Figure 14, indicate heavy 
long-wave activity at Santa Cruz and only light to moderate wave 
action at Monterey. The winds in Monterey Bay were northwest at 16 kts 
in the early morning decreasing to 2 kts in the evening. The surface 
air temperature averaged 55°F and clear skies were observed during the 
day. Although in contrast to the situation observed on 3 December 1969 
no relation of wind to long-wave activity can be drawn during this 
period. The 27.7 minute shelf-wave is readily apparent in both records 
with the 13.6 minute harmonic present at Santa Cruz alone. The Santa 
Cruz picture is quite similar to that experienced earlier on 3 December 
1969 at Santa Cruz. Monterey is generally weak in energy density and 
the shelf-wave appears to be the only active wave of any significance. 

14 January 1970 

The analyses of wave records on 14 January 1970 are the most dis- 
tinctive of the study. The sea surface traces, Figure 16, show moderate 
to heavy wave activity at both stations. During the period analyzed, 
winds were southerly at 35 kts and heavy rain was falling in the bay. 
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Sea Surface Heights (feet) Sea Surface Heights (feet) 





Figure 14 

Sea Surface Heights, Santa Cruz/Monterey 
2348 - 1848, 5-6 Dec 69 
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Frequency (millihertz ) 

Santa Cruz, 2348 - 1848, 5-6 Dec 69 




Frequency (millihertz ) 
Monterey, 2348 - 1848, 5-6 Dec 69 



Figure 15 

Spectral Wave Analysis, 6 Dec 69 
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Sea Surface Heights (feet) Sea Surface Heights (feet) 





Figure 16 

Sea Surface Heights, Santa Cruz/Monterey 
0810 - 0410, 14 Jan 70 - 15 Jan 70 



51 




Frequency (mi 1 1 iher t z ) 

Santa Cruz, 0810 - 0410, 14 Jan 70 - 15 Jan 70 




Frequency (millihertz) 

Monterey, 0810 - 0410, 14 Jan 70 - 15 Jan 70 



Figure 17 

Spectral Wave Analysis, 14 Jan 70 
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Winds decreased to 2 kts in the late evening, however, surface air 
temperature remained near 50°F during the entire period. The 27.7 
minute wave is the primary long-wave present at both locations. In 
addition the 18.9 minute wave is present at Santa Cruz but extremely 
weak at Monterey. Again the relationship between weather conditions 
and long-waves present is not characterized by any definite pattern. 

The 51.8 minute wave is present at both locations, but is weak. 

20 January 1970 

Wave activity at Santa Cruz was again more intense than that at 
Monterey. Winds were light and variable in the morning increasing to 
southerly at 15 kts with heavy rain in the early evening. The wave 
activity was similar to that experienced on 3 December 1969 showing 
moderate southerly winds, heavy wave activity at Santa Cruz and light 
activity at Monterey. However, with only two examples of this situ- 
ation, a clear pattern is not observed. The Santa Cruz spectrum shows 

some long -wave activity present at many discrete frequencies, however, 

* 

the 51.8 and 13.6 minute periods are the only recurring wave periods 
present. In this record the 27.7 minute shelf-wave was weak at Santa 
Cruz. The 36.1 and 27.7 minute waves were observed at Monterey, how- 
ever, overall wave activity at Monterey was weak. 

Looking at the overall picture, including Robinson's [1969] data, 
it was seen that the 27.7 minute wave evaluated as a shelf-wave, was 
the dominant long-wave in the bay. This 27.7 minute period is averaged 
over 5 frequency band widths therefore appears within the range of 26.3 
to 29.4 minutes. Various values of average depth of the shelf and 
distance to the shelf, in general, yield periods which fall within this 
range. The fundamental transverse period of oscillation was observed 
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Sea Surface Heights (feet) Sea Surface Heights (feet) 





Figure 18 

Sea Surface Heights, Santa Cruz/Monterey 
0215 - 2215, 20 Jan 70 
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Energy Density (meters 2 -sec) Energy Density (meters -sec) 

OOC 0G1 002 000 00 L 002 




Frequency (millihertz) 

Santa Cruz, 0215 - 2215, 20 Jan 70 




Figure 19 

Spectral Wave Analysis, 20 Jan 70 
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regularly at Monterey, but only once at Santa Cruz . The fundamental 
longitudinal period was significant at both locations most of the time. 
No definite relationship was drawn between atmospheric conditions and 
long-wave activity, however, a strong southerly wind seemed to lead to 
heavy wave activity in Santa Cruz. Also it was noted that most of the 
time, long-wave activity when present in the bay, was much more intense 
at Santa Cruz than at Monterey. 

C. INTERPRETATION OF CROSS SPECTRA AND PHASE 

Table VIII gives a summary of selected recurring long-period waves 
which have significant cross power when compared with other peaks in 
the same record. Some high energy-density long-waves which only appear 
once or twice in the entire study are excluded from the analyses. The 
cross spectra is a measure of the combined energy-density resulting 
from contributions of each power spectra but gives no information as to 
the forcing function of the wave or to whether the waves at the two 
locations are correlated with each other. Two isolated waves of the 
same period but operating independently, each having considerable 
energy, would yield a high cross power. Haase relationships determine 
whether or not the waves appearing at two locations are coupled. 

Figures 20 through 24 show the cross spectrum and phase differences 
of each long-wave between Monterey and Santa Cruz. Both raw spectra 
were smoothed over five frequency band widths. However, considerable 
difficulties arose when attempting to smooth the phase relationships. 

The process of smoothing over band widths assumes that the data is 
the result of a stationary random process. The long-waves analyzed in 
this study may be viewed as deterministic functions and therefore what 
is considered in the record is a number of deterministic functions 
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superimposed on so called n white noise". The ideal way to smooth the 
phase information would be to smooth over ensembles, however, this could 
not be done because the long -waves are, in general, a transient phenomenon 
and usually not a stationary process. The Fast Fourier Transform (FFT) 
method of calculating the power spectrum is relatively new and, while 
the raw phase information is considered correct, the correct method to 
smooth relationships when using the FFT have not been completely worked 
out [Enochson and Otmes , 1968]. 

In this study the raw phase information was calculated utilizing 
the cross spectra as a weight function such that 



e = 




i = 1,2, 3, 4, 5 



where 



e = smooth phase difference 
0 ^ = raw phase 

= raw cross power. 

This method was adopted after investigating several alternatives. The 
calculation resulted in those raw phase estimates associated with a high 
cross power to be dominant in the smoothing process. Those relation- 
ships are shown in Figures 20 through 24 and summarized in Table IX for 
selected periods. There are reasons to doubt the accuracy of calcu- 
lating the smoothed phase using FFT in this manner which circumvents 
computing the correlation function. Hence, the smoothed phase calcu- 
lation must be viewed with reservations. 

Table VIII is a summary of significant cross spectral peaks for 
selected long-waves which recur throughout the five days analyzed. 
Figures 20 through 24 show the raw and smoothed cross spectra for the 
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PERIOD 

(minutes) 


6 Nov 69 


3 Dec 69 


DATE 

6 Dec 69 


14 Jan 70 


20 Jan 70 


51.8 


X 




X 


X 




42.5 


X 










36.1 




X 






X 


27.7 






X 


X 


X 


22.5 












18.9 






X 


X 




16.3 


X 










13.6 




X 








10.1 


X 


X 


X 






9.3 






X 










TABLE 


VIII 







Summary of Significant Cross Spectral Peaks in Monterey Bay 
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Figure 20 

Cross Spectra and Phase Relationships, Monterey/Santa Cruz, 6 Nov 69 
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Cross Spectra and Phase Relationships, Monterey/Santa Cruz, 3 Dec 69 
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Cross Spectra and Phase Relationships, Monterey/Santa Cruz, 6 Dec 69 
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Cross Spectra and Phase Relationships, Monterey/Santa Cruz 

14 Jan 70 
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Cross Spectra and Phase Relationships, Monterey/Santa Cruz, 20 Jan 70 



dates analyzed. As expected the cross spectra is dominated by the 27.7 
minute shelf-wave and its associated harmonics. The cross spectral 
summary is in general agreement with Table VII, the summary of indivi- 

i 

dual spectral peaks. The phase differences, between a particular wave 
at Santa Cruz and Monterey appear to be random, which implies the 
Monterey Canyon acts as a barrier dividing the bay into separate north- 
south basins which oscillate independently. This is in general agreement 
with Wilson [1965]. 

D. COMPARISON WITH WILSON'S THREE-DIMENSIONAL NUMERICAL MODEL 

Wilson [1965] assumed a boundary nodal line for his solution of his 
three-dimensional model of the oscillating characteristics of Monterey 
Bay. This node is a line drawn from Pt . Santa Cruz to Pt . Pinos on the 
Monterey Pens inula. This assumption tends to be the greatest objection 
to Wilson's results, although an enclosed bay will usually have a node 
across the seaward opening for particular modes of oscillation. The 
oscillating characteristics are depicted in Figures 25 through 28 where 
Figure 25(a) gives the general bathemetry of the bay. These depths 
were used as grid points in order to generate the numerical solution. 

The other figures show the increasing complex longitudinal modes of 
oscillation. The contour lines are the amplitudes of the water level 
normalized to the highest anti-node for the mode. Each complex mode 
of oscillation is drawn in a simplified version to the left of the 
figure. The nodal line assumption is particularly constraining to the 
longitudinal mode of oscillation which can be seen in Figure 25 where 
mode 1 (i.e., the fundamental) tends to oscillate transversely due to 
the nodal constraint. The lowest modes would be most affected by the 
nodal assumption and would have less constraint on higher modes. 
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Numerical Calculations of Modes of Oscillation of Monterey Bay (from Wilson, 1965) 
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Figure 26 

Numerical Calculations of Modes of Oscillation of Monterey Bay (from Wilson, 1965) 
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Numerical Calculations of Modes of Oscillation of Monterey Bay (from Wilson, 1965) 
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Numerical Calculations of Modes of Oscillation of Monterey Bay (from Wilson, 1965) 



* 











The periods of oscillation of the longitudinal mode were calculated 

as, 

! T n = 44 ’ 2 ’ 29 - 6 > 28 - 2 > 23.3, 21.6, 20.4, 

19.4, 18.7, 17.6, mins. 

Mode 1 has a fundamental period of 44.2 minutes and indicates strong 
oscillation at Santa Cruz and weak at Monterey. Mode 2 is similar to 
mode 1 while mode 3 gives the only strong oscillation at Monterey and 
weak oscillations at Santa Cruz. Mode 4 is similar to modes 1 and 2 . 
but is somewhat more complex as are successive modes. 

Wilson does state that his assumed node may very well be incorrect 
and possible should be moved further seaward. The boundary conditions 
of the model would then be changed, resulting in new oscillating 

periods and characteristics. This would result in a lengthening of 
the fundamental period. 

Wilson also computed two-dimensional modes of oscillation of 

Monterey Bay. This is the transverse oscillation of the east-west 

extremeties bounded at Moss landing on the east and the assumed node on 

the west. The sequence of modal periods of oscillation were calculated 
as 

T „ = 32.3, 14.3, 9.5, 7.0, ... mins. 

The profiles of water-surface elevation for the transverse mode are 
shown in Figure 29. 

Considering both the longitudinal and transverse modes of oscillation, 
Wilson notes that decreased solution confidence is placed on each suc- 
cessive harmonic, however, the highest mode in which he has confidence 
is not stated. This increasing error in successive modal solutions can 
be decreased by choosing a finer grid. Wilson's model probably is most 
valid for modes greater than mode 2 and less than mode 7. 
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As stated earlier, Wilson concluded that the deep Monterey Canyon 
does exert a profound influence on Monterey Bay oscillations as it 
effeptively sets a barrier to free oscillations, between the north and 

j 

south portions of the bay. The result is that wave activity acting in 
both locations would remain largely uncoupled and would give rise to a 
random phase. 

Table IX summarizes some of the results of this study with that of 
Wilson's model. It is apparent that the long-wave activity at Santa 
Cruz is more intense than that at Monterey. This tends to support the 
conclusion that the two basins act independently. Robinson [1965] 
found similar magnitudes at both stations and concluded that the 
oscillations were not uncoupled. The fundamental longitudinal mode of 
51.8 minutes was found to be weak at Monterey and stronger at Santa 
Cruz thus supporting Wilson's findings. The 20.5 minute wave corres- 
ponding to Wilson's mode 6 is in general weak at both stations. This 
is in agreement with Wilson's results. The greatest contrast between 
the results of this study and that of Wilson's is the effect and identi- 
fication of the 27.3 minute wave, evaluated as a shelf-wave here, but 
as the third mode of oscillation by Wilson. This wave, according to 
Wilson, should be strong at Monterey and weak at Santa Cruz. In general, 
this was not found to be the case. The ratios of energy-density at 
Santa Cruz to energy-density at Monterey is greater than 1 on 6 November, 
3 December and 6 December indicating the opposite of Wilson's findings. 
The last two periods studied show a ratio of .94 on 14 January 1970 and 
.49 on 20 January 1970, neither of which support Wilson's calculated 
results. It is considered that the 27.3 minute wave predicted by Wilson 
is in fact present during long-wave oscillations, but is evaluated as a 
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shelf-wave vice a harmonic of the fundamental frequency of oscillation 
of the bay. Most of the other periods predicted by Wilson were found to 

exist. However, it must be remembered that a discrete Fourier analysis 

| 

was performed and the data smoothed. This has the result of "forcing 1 * 
some energy into discrete band widths allowing only certain periods to 
appear in the spectrum. 

As a result of this comparison the present study tends to verify 
many of Wilson's conclusions. The observations that wave amplitudes 
are generally higher at Santa Cruz and that the phase relationships tend 
to be random lead to the conclusion that oscillations of Monterey Bay 
at Santa Cruz and Monterey tend to operate independently of each other 
due to the impedence barrier formed by Monterey Canyon. 



74 



V. ERROR ANALYSIS 



Errors encountered during the course of the study were basically 
of three types: 

(1) Time errors can be introduced when recording raw data on an 
analog trace. 

(2) Errors may be generated in matching the Santa Cruz raw 
data to Monterey data. 

(3) The effect of the numerical methods of calculating the 
power-spectrum and the cross spectrum can lead to mis- 
interpretation of results. 

Each of these effects are examined below in order to determine, in a 
qualitative sense, their effect on the calculated results. 

Time Induced Errors 

The tide gages at Santa Cruz and Monterey are both mechanical 
devices which must be wound at required intervals. Both were maintained 
within the required periods. The sampling interval of 34.495 seconds 
for Monterey and 6.00 seconds for Santa Cruz were computed by inches/ 
time relationships between two successive time checks on each record 
for each period analyzed. The times at Santa Cruz were recorded from 
radio station WWV (Naval Observatory Time) and are correct to the second. 
The times at Monterey were obtained from the telephone company and are 
accurate within + 5 seconds. Time checks were maintained daily for the 
entire period. The recording speed at Monterey was found to be 1.03 
+ .01 inches/hour and 6.00+ .02 inches/hour for Santa Cruz. Synchronized 
records are critical to phase computations. However, the resulting 
error in obtaining the raw data was considered insignificant. 
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Errors Generated in Aligning Records 



The Santa Cruz record was aligned with the Monterey record by 
converting the coordinate system from curvilinear coordinates to 
rectilinear coordinates. Linear interpolation was then used to match 
data points. 

The Santa Cruz data yielded a data point every six seconds giving 
high resolution in digitizing. The original 12,000 plus data points 
were reduced to 2048 points. Reduction and conversion of data in this 
manner resulted in very minor raw data errors. 

Numerical Calculation Methods 

The effect of approximating a continuous tidal trace with a step- 
function and transforming the step-function by numerical methods through 
use of a digital computer has several effects on the final analyses, two 
of which could result in misleading results of interpretation. 

When a continuous curve is approximated by a step function, a 
secondary wave form of high frequency energy is superimposed on the 
recorded wave form. To evaluate this possible error the sampling rate 
was reduced and records were analyzed and checked for aliasing. It was 
found the aliased power was quite small and in no case greater than 10$. 

The second source of error is the effect of transforming the record 
into discrete values of energy-density as a function of discrete periods. 
Thus, the energy is forced into certain band widths and centered about 
the middle frequency. As a result of smoothing over five frequencies, 
the frequency interval was .00007 cps ranging in periods from 397.4 
minutes to 1.16 minutes. Comparison of raw and smooth power and cross 
spectra show that the raw spectra was not smeared due to smoothing the 
original 1025 spectral estimates as "spikes 11 in the raw spectra were 
still present after the smoothing process was completed. 
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VI. SUMMARY 



i 

Power spectra analyses and cross spectral analyses of wave records 
in Monterey Bay indicate the presence of persistent unique periods 
with the phase differences between Monterey and Santa Cruz tending to 
be random. Coupled with this and the fact that wave energy at Santa 
Cruz is in general greater than that at Monterey, it is concluded that 
the Monterey Submarine Canyon has a profound effect on seiching motions 
within the bay. This effect appears to divide the Monterey Bay into 
north and south basins where oscillations occur essentially independent 
of each other. 
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APPENDIX A 



DIGITIZING GRAPHICAL ANALOG RECORDS USING THE CAIMA COMPANY MODEL 480 

The Calma digitizer reduces analog graphical data to digital 
magnetic tape for computer processing and analysis. It consists of a 
freestanding tracing bed and a separate electronics/recorder module. 

The tracing bed is motor-driven for adjustment to the most comfortable 
digitizing position . 

To digitize analog graphical data directly on magnetic tape, the 
operator traces the analog curve with a moveable stylus/carriage assembly. 
As the curve is traced, movements of the stylus in either the X or Y 
direction causes pulsed to be generated by the direct digital pickoff 
mechanism. Incremental motions are recorded as characters on magnetic 
tape, leaving the task of summing the increments for whole-value coordi- 
nates to the final processing computer. This method of digitizing is 
both faster and more accurate than template methods. 

The digitizer reads and records data in the X and Y direction, with 
a sampling interval which can be set to 0.01, 0.02, 0.04, 0.08, or 0.15 
inches. The maximum absolute sampling error for the machine is 0.012 
inches. The output is external BCD, stored on 556 bpi, 7 channel, tape. 
The tape can be made compatible with the IBM/360 system, however, it is 
easier to use the CDC/6500 computer. 

A FORTRAN IV program wirtten for the CDC/6500 can be utilized for 
interpreting the digitized record. 

The program accomplished the following: 

(1) Converts the data from external BCD to display code. 

(2) Arranges the data in column matrices of 80 character length. 
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(3) Interrogated each character in the matrix to determine if 
the character is a 

(a) Point flag - a designator to label a specific data 
point of interest 

(b) Plus or minus travel in the X-direction. 

(c) Plus or minus travel in the Y-direction. 

(4) Sums the Y direction travel. 

(5) Records incremental travel in X direction (plus or minus). 
This allows for negative travel while digitizing. Negative travel 

is common if digitizing a curvilinear trace of if variations in data are 
closely spaced. 

(6) Prints the Y value V(N), for each designated increment in 
the X direction. 

(7) Prints the X increment value U(N), showing the incremental 
travel (plus or minus) for which the Y value was obtained. 

(8) Punch cards for V(N) and U(N) which are used as inputs for 
other computer programs as desired. 

Headings may be entered on the tape with keyboard control, but 
CONVERT is not designed to read tape headings nor does the program have 
the capability to search for a particular set of data on the tape. Con- 
sequently, all data on the tape will be analyzed each time CONVERT is 
used. 

The following sequence of operations can be used as a guide for the 
digitizing procedure: 

(1) Mount chart on tracing bed so that horizontal axis of chart 
is aligned with X axis of digitizer. 

(2) Turn system power on. 
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(3) Load tape. 

(4) Initiate "Load Forward" on tape module. 

(5) Press RECORD ERROR. (If recorder is properly loaded, RECORD 

i 

ERROR indicator will go out when push button is pressed,... 
digitizer will automatically shift to KEYBOARD MODE.) 

(6) Enter required identifiers from keyboard if necessary. 

(Not required with CONVERT program.) 

(7) Position stylus at beginning of curve. Lock Y-axis, and 
check alignment. 

(8) Initiate TRACER Mode, de-activate SKIP (SKIP button out, 
light off), and unlock stylus. 

(9) FIAG the tape once. 

/ 

(10) Trace the record as desired. The first point which will not 
be computed, or printed, will be X = 0, Y = 0. The speed at 
which the stylus can be moved without destroying data 
accumulated is well in excess of the speed normally used to 
digitize the record. However, if exceeded, the RECORD ERROR 
will light and record will be deleted. 

(11) When tracing is complete, lock the stylus, initiate FLAG 
PEDAL twice and IRG button once. 

(12) Press rewind button on tape module and remove tape. 

The user of the Calma digitizer can reap the benefits of the instru- 
ment’s speed and accuracy only if the digitizing operator understands 
and applies the system’s basic rules of operation. It is strongly 
suggested that new digitizing operators read, study and understand the 
textual instructions (CA1MA MODEL 480 DIGITIZER, Instruction and Main- 
tenance Manual) before attempting to operate the digitizer. 
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(13) The magnetic tape from the digitizing process is used as INPUT 
to program CONVERT. The approximate number of words on the 
tape should be estimated so that CONVERT may be dimensioned 
accordingly . 

In addition, if the tape contains any parity errors, or was not 
deguassed properly prior to usage, the program will only interpret a 
portion of the data or possibly not run at all. These problems have 
been encountered in the past. 

(14) The program CONVERT is included and is self-explanatory. 

Output consists of (a) sampling rate, (b) length of record 
(sec), (c) amplitude values (inch) V(N), and (d) incremental 
travel in X-direction U(N). These may be summed in another 
program to determine proper time or spatial frame. 

(15) It must be remembered that CONVERT is FORTRAN IV written for 
the CDC/6500 located at FNWC Monterey, and must be adapted 
for other computers. 
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PROGRAM CONVERT! I NPUT , OUTPUT , PUNCH ) 

THIS PROGRAM CONVERTS THE DIGITIZED RECORD OF THE 
MONTEREY TIDE GUAGE FROM EXTERNAL BCD TO AMPLITUDE 
VALUES FOR A CONSTANT SAMPLING PATE DELT, AVGX IS THE 
AVERAGE DISTANCE THE RECORDING CRUM ADVANCES PER ONE 
HOUR PEAL TIME, I, E, THE DATA SAMPLING RATE. 
REFERENCE-CALMA CO, MODEL 303 DIGITIZER INSTRUCTION 
MANUAL. DIMENSION STATEMENT MAY BE MODIFIED TO 
COMPUTE AS MANY DATA POINTS AS REQUIRED. 

DIMENSION U( 17000) ,V( 17000) ,N( 80 ) , I BUFF( 17000) ,NK(80) 
PROGRAM READS L SETS OF DATA FROM TAPE 
READ 97, L 

97 FORMAT (12) 

DO 116 J J= 1 » L 
READ 98, AVGX 

98 FORMAT! F10.3) 

COMPUTE SAMPLING RATE (SECONDS) 

DELT=36.0/AVGX 
PRINT 99, DELT 

99 FORVAT( 1H0,2X,25H SAMPLING INTERVAL EQUALS , IX , E 15 .7 , IX 
1 ,7H$EC0NDS ) 

PRINT 700, JJ 

700 FORMAT! 1H0,9HJ J EQUALS, 2X,I2) 

ZERO FILL AMPLITUDE ARRAYS 
DO 100 1=1,17000 
U( I ) = 0 . 0 
V( I ) = C.C 

100 CONTINUE 

ZEROIZE INPUT/CUTPUT BUFFER 

101 DO 202 1=1,17000 
I BUF F ( I ) =0 

202 CONTINUE 
COUNTX=C . 0 
COUN TY = 0 . 0 
NUM= 17000 

COUNT NUMBER OF DATA POINTS ACCUMULATED 

M= 1 

K8=0 

CALL SUBROUTINE LIOF TO READ INPUT TAPE WHERE NPAR IS 
A PARITY CHECK, AND NEOF CHECKS FOR END OF FILE 
CALL L I OF ( 5LRBCD1, I BUFF, NUM, NPAR, NEOF) 

IF(NEOF) 602,200 

200 K=- 7 
KF = 0 

201 K=K+ 8 
KA = K + 8 
KC = 0 

DO 102 KB=K,KA 

irr — is r x i 

NK ( KC ) = I BUFF ( KB ) 

IFINK(KC).EO.O) KF=l 

102 CONTINUE 

DECODE TAPE SO WORDS AT A TIME 
DECODE! 80, 103, NK) ( N ( I ) , I = 1 , 80 ) 

103 FORMAT ( 80R1 ) 

INTERPRET BCD ON TAPE AND COMPUTE AMPLITUDES 

DO 104 19=1,80 

IF(N( 19). EG. SOB) GC TO 106 

I F ( N ( I 9 ) . EQ. 55B ) GO TO 107 

IF! N ( 19) .EQ.34B) GO TO 103 

I F ( N ( 19) .EQ.478) GO TO 105 

SYMBOL / ( 50B ) REPRESENTS AN INCREMENT TRAVEL IN THE 
MINUS X OR Y DIRECTION BY THE DIGITIZER. 

SYMBOL 0, ( 55B) , REPRESENTS A ZERO INCREMENT TRAVEL IN 
THE X OP Y DIRECTION BY THE DIGITIZER. 

SYMBOL 1,( 34B) , REPRESENTS AN INCREMENT TRAVEL IN THE 
POSITIVE X OR Y DIRECTION BY THE DIGITIZER. 

SYMBOL *,(47B), IS A FLAG INSERTED IN THE RECORD BY 
THE PERSON DIGITIZING. 

GO TO 104 

105 PRINT 205, M, 19 
205 FORMAT! 1H0, 21 10) 
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C THIS SECTION ALLOWS FIRST TEN WORDS OF TAPE TO BE 

C IDENTIFICATION DATA IF OESIRED. 

IF(M.GT.IO) GO TO 113 
GO TO 104 

C COMPUTE NEGATIVE TRAVEL 

l:06 R X=- 0 . C 1 
K8=K8+ 1 
GO TO 109 

C COMPUTE ZERO TRAVEL 

107 RX=0 .0 
K6=K 8 + 1 
GO TO 1C9 

C COMPUTE POSITIVE TRAVEL 

108 RX=0 . 0 1 
K8 = K 8 + 1 

109 K3 = K 8/2 
K3=2 *K 3 

IF(K3.E0.K8) GO TO ill 
C01JNTX=CCLNTX+RX 
I F { CCUNTX.NE .0.0) GO TC 110 
GO TO 104 

C COMPUTE X INCREMENT TRAVEL 

110 U( M)=CCLNTX 
GO TO 104 

111 COUNTY =CCUNTY + RX 
IF(COUNTX.NE.O.O) GO TC 112 
GO TO 104 

C COMPUTE AMPLITUDE VALUE 

112 V ( M ) =C CLNT Y 
C0UNTX=0. 0 
M=M+ 1 

C STOP IF DATA POINTS EXCEED ARRAY 

IF( M.GT. 17000) GO TO 600 

104 CONTINUE 

IF(KF.EQ.l) GO TO 113 
GO TO 201 

C TOTAL TIME OF THE RECORD EQUALS THE NUMBER OF DATA 

C POINTS TIMES THE SAMPLING I NTE RV AL , DELT . 

T I M E = ( M*DELT)/3600.0 
PRINT 115, TIME 

115 FORMAT ( 1H ,10X,28H TOTAL TIME OF RECORD EQUALS , 2X , E 1 5 . 

1 7, 2X, 7H HOURS,) 

C PRINT AND PUNCH AMPLITUDE VALUES 

PRINT 215, ( V( I ) , I =1 , M) 

PUNCH 213, (V( I ) ,1 = 1, M) 

C PRINT AND PUNCH X INCREMENTS 

PRINT 215, <U( I ) ,I=1,M) 

PUNCH 215, (U( I) ,1=1, M) 

215 FORMAT ( 1H ,10X,14F7.2) 

213 FORMAT ( 14F5.2) 

116 CONTINUE 
STOP 

600 PRINT 601 

601 FORMAT ( 1H1 »20X,35H**** U AND V SPACE INADEQUATE ***♦*) 

602 STOP 
END 
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APPENDIX B 



COMPUTATION OF FAST FOURIER TRANSFORM USING IBM/360 SUBROUTINE RHARM 
This subroutine is a one-dimens ional Fast Fourier Transform (FFT) 
analysis based on the Cooley and Tukey [1965] algorithm. The program 
is designed to analyze I data points where 

I = 2 M H « 3, 4, 5, . . . , 20. 

The FFT not only greatly reduces the number of calculations from 
earlier analysis schemes, but also reduces the round-off errors in the 
Fourier coefficients. The computations and round-off errors are 
reduced by log 2 (I)/l- 

The subroutine HARM is then called by RHARM to compute the complex 
coefficients : 

■ 5 £ < X 2j- IX 2 j+ l> ; 

Then for K = 1,2,...,| - 1, | (with the 
A K = 2 + » A K = 

and , 

C K - I W K + ^ ' » 

C N-K - I - A K ^ 

for K = 1,2,. . . - 1. 



N 






K = 0,1,2, ...N-l. 

bar denoting conjugation) 



2 ^N-K"^ 



K). 



, i, _ i o N i N 

for K - 1,2 , . . «2 “ 1, 2 



- ~) 
20 
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If we let Cq = ^ (ReA^ - ImA^) 

C N * 1 (ReA 0 + In V 

We finally have: 

A q = 2 Re (C 0 ), B k = -21m (C K > , 

B 0 = °> Ajj = 2 Re (C N > , 

A K = 2 Re (C R ), B n = 0. 

The Fourier coefficients = 0, A., B.....A.., B„ = 0 are 

0 0 l i N N 

obtained for input X ^ , j = 0, 1 ,2 , • . .2N-1 , where N = 1/2, for the 
following equation. 

N-l 

Xj = ^ A Q + E (A^ COS (ttJK/N) + B R SIN (TTjk/N) + 

K= 1 

5 *« 
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APPENDIX C 



I 

LEAST SQUARES CURVE FITTING 

The scheme used to accomplish the least squares fit for a sixth- 
degree polynomial is described below. Given a set of data pairs (X^ , Y^) 
(i = 0,N) which can be interpreted as measured values of the coordinates 
of the points on the graph of y = f(x), assume that the unknown 
function f(x) can be approximated by a linear combination of suitably 
chosen functions, fQ(x), f^(x), ... f^(x) t * ie ^ orm **(x) = A Qfg( x ) + 
A^f^(x) + A 2 f 2 (x) + A^f^Cx) + A^f^ (x) + A^f^(x) + A^f^(x) where the 
unknown coefficients A^, A^, . .., A^ are independent parameters to be 
determined and the degree of least squares polynomial, M is such that 
M = 6 < N. The difference between the approximating function value 
F(xj) and the corresponding data value Y ^ , is called the residual r 
and is defined by the relation 

rj = F(X.) - Y. (i = 0, N). 

The function F(x) that best approximates the given set of data in 

a least squares sense is that linear combination A^ fg(x) + + 

... + A.f^Cx) of functions f. (x) that produces the minimum value of 
oo k. 

the sum Q of the squared residuals where 

Q = S r. 2 h S [f(X.) - Y. 2 ]. 
i i 

For the case of curve -smoothing by polynomial least squares, we 
approximate the function y = f(x) over the range of data points 
(X^, ^i^ ( i = N). The parameters Aq, A^, . . . A^ are then deter- 

mined such that 
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« ' E r i 2 * E [W - Y i 2 ] 

1 1 

is £ minimum. That is, a 6th degree polynomial curve is fitted to the 

i 

data points in a least squares sense as defined earlier. The normal 
equations for the least-squares polynomial can be written as 
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where ; 

N =* number of data pairs 
= time increment (seconds) 

= amplitude values (inches) 

The solution of the matrix is the sixth-degree polynomial, 

V*> = A q + A^X 1 j+ A 2 X 2 2 + A 3 X 3 3 + A 4 X 4 4 + A 5 X 5 5 + A 6 X 6 6 

which is then subtracted from each data point to remove tidal effects, 
acting as a high-pass filter. Choosing a sample interval of 90 
increments, 23 data pairs were utilized for the least squares fit. 
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